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Abstract 

The aim of this paper is to understand the tendency to organization of the 
turbulence in two-dimensional ideal fluids. We show that nonlinear processes 
as inverse cascade of the energy and vorticity concentration are essentially de¬ 
termined by trajectory trapping or eddying. The statistics of the trajectories 
of the vorticity elements is studied using a semianalytic method. The sepa¬ 
ration of the positive and negative vorticity is due to the attraction produced 
by a large scale vortex on the small scale vortices of the same sign. More pre¬ 
cisely, a large scale velocity is shown to determine average transverse drifts, 
which have opposite orientations for positive and negative vorticity. They 
appear in the presence of trapping and lead to energy flow to large scales 
due to the increase of the circulation of the large vortex. Recent results on 
drift turbulence evolution in magnetically confined plasmas are discussed in 
order to underline the idea that there is a link between the inverse cascade 
and trajectory trapping. The physical mechanisms are different in fluids and 
plasmas due to the different types of nonlinearities of the two systems, but 
trajectory trapping has the main role in both cases. 

Keywords: test particle statistics, turbulence, self-organization, vorticity, 
Euler fluid 


1. Introduction 

The two-dimensional fluid turbulence represented during more than 70 
years an active field of research (see the review papers jl|, 0], (3[ and the ref¬ 
erences therein). It has many applications in different areas as fluid dynam¬ 
ics, meteorology, oceanography, fusion plasmas, superfluids, superconductors 
and astrophysics, in spite of the fact that it provides only idealized models 
for physical systems that are always three-dimensional. 
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The two-dimensional turbulence has a self-organizing character, which is 
related to the invariance of both energy and enstrophy in ideal (inviscid) 
fluids. Numerical studies of the decaying turbulence clearly show this prop¬ 
erty and the associated scaling behaviour i, 0, 0. The enstrophy has a 
direct cascade (to small scales), but with a complex evolution characterized 
by the presence of inverse cascade in isolated regions 0, 0. The energy 
has an inverse cascade (to large scales) that leads to the emergence of large 
quasi coherent vortices. The process of self-organization can continue until 
the coherent vortices reach the size of the system 0, [13 • This behaviour 
was explained in the representation of point-like vortices by a negative tem¬ 
perature (0, 0. 0) or by the property of self-duality of the associated 
field theoretical model (jl3j|). The latter approach was extended to a model 
of planetary atmosphere and magnetized plasmas [l3 |. 

This paper deals with the turbulent states and studies the self-organization 
during its initial stage, before the emergence of large coherent vortices of the 
system size. 

We show that trajectory trapping or eddying in the structure of the tur¬ 
bulence is the main physical reason for the strong nonlinear effects that are 
observed in two-dimensional ideal fluids. This conclusion is drawn from a 
study of the statistics of test particles (tracers) in turbulent Euler fluids. 

This study is based on a series of recent results on the statistical prop¬ 
erties of test particle trajectories in incompressible two-dimensional velocity 
holds. Numerical simulations have shown that trajectories are complex, as 
they have both random and quasi-coherent aspects. A typical trajectory, 
shown in Figure 1, is a random sequence of long jumps and trapping events 
that consists of winding on almost closed paths. Analytical methods that 
describe the statistics of these trajectories were developed [13, [3] and used 
for understanding various aspects of turbulent transport. It was shown that 
they provide a very good description of the nonlinear effects produced by 
trajectory trapping or eddying and reasonably accurate quantitative results 
for the diffusion coefficients and for other statistical averages. 

The conclusion of these studies is that trajectory trapping or eddying 
leads to nonstandard statistics: memory effects (represented by long time 
Lagrangian correlation), strongly modified transport coefficients and non- 
Gaussian distributions of displacements. It was also shown that trapping 
determines a large degree of coherence in the sense that bundles of trajec¬ 
tories that start from neighboring points remain close for very long time 
compared to the eddying time. Trapped trajectories form quasi-coherent 
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structures similar to fluid vortices. Extensive theoretical 17]- 2l| and nu¬ 
merical studies ( 22 }] -[HI have contributed significantly in the last decades to 
the understanding of the turbulent transport in laboratory or space plasmas, 
in fluids or in stochastic magnetic fields. 

The connection between the test particle trapping and the turbulence 
evolution was discussed in 25] based on a study of test modes on turbulent 


plasmas. Analytical results that are in agreement with numerical simula¬ 
tions were obtained, and they allowed to deduce a new physical perspective 
on the nonlinear process of generation of large scale correlations (inverse cas¬ 
cade) and of zonal flow modes. Essentially, they are effects of ion trajectory 
trapping or eddying. 

We show here that trapping has an essential role in two-dimensional fluid 
turbulence. A nonlinear effect produced by trapping of the vorticity elements 
explains the separation of positive and negative vorticity and the inverse 
cascade of the energy. 

The paper is organized as follows. The problem of test particle or tracer 
transport is defined in Section 2.1 and Section 2.2 contains a short pre¬ 
sentation of the analytical statistical approach, the decorrelation trajectory 
method (DTM). The effects of trapping or eddying on tracer transport and 
on the statistical characteristics of the trajectories are discussed in Section 3. 
This section contains a review of the previous work and new results on the 
modifications of the trajectory structures determined by an average velocity. 
The effects of trapping on the decaying two-dimensional turbulence in ideal 
fluids are analyzed in Section 4.1. The physical process that determines the 
separation of the positive and negative vorticity and leads to the inverse cas¬ 
cade of the energy is identified. The average speed of vorticity separation 
is estimated. Section 4.2 is a short discussion on recent results 011 plasma 
turbulence evolution, which show that trajectory trapping has an essential 
role in the inverse cascade, although the physical mechanism is completely 
different. The conclusions are summarized in Section 5. 


2. Test particle transport and the statistical method 

2.1. The problem 

The problem of test particle or tracer advection in two-dimensional in¬ 
compressible velocity fields is described by the stochastic equation: 

dx(f) 


dt 


= v [x(t),i] + V d e y , 


( 1 ) 
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where x(t) represents the trajectory in the plane (e^, e y ) and V^ey is an 
average velocity that is taken in the e y direction. The stochastic velocity 
v(x, t ) is incompressible [V ■ v(x, t ) = 0] and is represented by a scalar held, 
the stochastic potential or stream function 

v(x,t) = e 2 x V0 = (~d y <l), d x <t>). (2) 

The potential </>(x, t) is considered to be a stationary and homogeneous Gaus¬ 
sian stochastic held, with zero average and given two-point Eulerian correla¬ 
tion function (EC) 


E(x, t ) = (0(x', t') Vj(x' + x, t' + t)) (3) 

where (...) denotes the statistical average over the realizations of </>(x, t) or 
the integral over x' and t' . The main parameters of the EC are: the amplitude 
of the potential huctuations /3 2 = £7(0, 0), the correlation length A c and the 
correlation time r c , which are the characteristic length and time of the decay 
of the function E(x, t ). The EC’s of the velocity components are obtained as 
space derivatives of E(x,t) and the amplitude of the stochastic velocity is 
= 0 /\ c . 

Starting from the above statistical description of the stochastic potential 
and from an explicit EC one has to determine the statistical properties of 
the trajectories. This problem is nonlinear due to the space dependence of 
the potential, which leads to x dependence of the EC (EJ)- 

The trajectories are solutions of a Hamiltonian system and thus, for time 
independent potential 0(x), their paths are the contour lines of the total 
potential c/>(x) + xVd (the stream lines). The velocity is a continuous function 
of x and t in each realization and it determines an unique trajectory as the 
solution Eq. (jT]) with the initial condition x(0) = 0 . 

For Vd = 0 the stream lines are closed curves (with the exception of 
one contour line, 0(x) = 0, which extends to infinite), and the trajectories 
are periodic functions of time. They correspond to permanent trapping. 
This invariance property is approximately maintained for potentials that are 
weakly perturbed by a slow time variation. The trajectories approximately 
follow the contour lines of 0(x, t) on almost closed paths for time intervals 
that are larger than the time of flight Tfi = \ C /V. This leads to typical 
trajectories similar to the example shown in Figure 1. They have a complex 
structure that consists of a random sequence of trapping events and long 
jumps. Trapping appears at different scales when the particles are close 
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to the maxima or minima of the potential. The large displacements are 
produced when the trajectories are at small absolute values of the potential. 
The importance of trapping is measured by the Kubo number defined by 



Figure 1: A typical trajectory obtained from Eq.(l) for A' = 10, Vd = 0 and the stream 
function used in the simulations presented in [30]. 


T, 


K = — = 

r fi 


Vr c 

A, 


(4) 


Trajectory trapping appears for K > 1 and becomes stronger as K increases 
up to the limit of static or frozen fields (. K , r c = oo) when the trapping is 
permanent. 

An average velocity Vd changes the topology of the total potential 0(x) + 
xVd by generating bunches of opened contour lines along its direction. When 
r = Vd/V < 1, islands of closed contour lines persist between the bunches of 
opened lines, which enables trajectory trapping. At large average velocities 
r 3> 1 all the contour lines are opened and test particle cannot be trapped. 
The average velocity defines a characteristic time for traversing the corre¬ 
lation length with the average velocity, r* = A C /Vd, and a dimensionless 
parameter similar to the Kubo number 


AT = 


T fi 


V_ 

Vd 


1 

r 


(5) 


The condition that trapped trajectories are obtained from Eq. © is 
K > 1 and AT >1 (or r < 1). 
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Test particle transport is essentially determined by the Lagrangian veloc¬ 
ity correlation (LVC) defined by 

Uj(t) = (vi [x(0), 0] vj [x(t), t ]} (6) 

for a stationary process. As shown by Taylor (26|, the mean square displace¬ 
ment (xf(t)) and its derivative, the running diffusion coefficient Di(t ), are 
integrals of the LVC 


{x 2 (t)) = 2 [dr Lair) (t - r), 

Jo 

(7) 

Di(t ) = [ dr Lh(t). 

Jo 

(8) 

The asymptotic values of D,(f) 


poo 

Xi = lim Di(t) = / Ljj(r)dr 

i-KX> J Q 

(9) 


are the diffusion coefficients. The time dependent diffusion coefficients pro¬ 
vide the ” microscopic” characteristics of the transport process. The diffusion 
at the transport space-time scales, which are much larger than X t and r c , is 
described by the asymptotic values x%- 

The scaling of the diffusion coefficient in the Kubo number can be ob¬ 
tained by a simple estimation based on the general shape of the LVC (!hj). It 
is usually a function that decays to zero from the value V 2 at t = 0 in a char¬ 
acteristic time Td (the decorrelation time). For small Kubo numbers K -C 1 
(t c <C Tfi), the time variation of the velocity field is fast and the particles 
cannot ’’see” the space structure of the velocity field. In this quasilinear 
regime (or the weak turbulence case) = r c and the diffusion coefficient 
® is X q i ~ V 2 t c = ( A 1/t c )K 2 . In the nonlinear regime K >■ 1 ( 77 / r c ), 
the decay time of the LVC is Td — Tfi and the diffusion scales as the Bohm 
diffusion coefficient \b ~ V’V/j — ( A 1/t c )K, which does not dependent on 
t c . The shape of the EC determines only a numerical factor in the diffusion 
coefficient. 

The Bohm scaling does not apply in the case of diffusion in 2-dimensional 
divergence-free velocity fields described by Eq. ©■ Numerical simulations 
have shown that the diffusion coefficient does not saturate as r c increases, 
but it decreases with the increase of the correlation time at K > 1 (r c > Tfi). 
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Trajectory trapping or eddying was found to produce this effect [27[. 
scaling in the nonlinear regime is 


Xtr » (>i/ T c)K“ = PK' 


a—1 


The 


( 10 ) 


where a < 1. In particular, in the limit of static or frozen potential K —> oo, 
the process is subdiffusive with Xtr 0- 

2.2. The decorrelation trajectory method (DTM) 

Several methods (as Corrsin approximation |28j], (jjL direct interaction 
approximation [30i], the renormalization group technique |3l|) were developed 
for determining the LVC corresponding to given EC. All these methods lead 
to finite diffusion coefficient of Bohm type in the limit K — y oo, which shows 
that they are not adequate for the two-dimensional incompressible velocity 
fields. This process was studied especially by means of direct numerical sim¬ 
ulations (|32J and the reference there in) or by developing simplified models 


33]. There is a theoretical estimation [34| based on the analogy with the 


fractal structure of the landscapes. Using the theory of percolation, it finds 
the scaling of the diffusion coefficient of the type fjTOjl with a = 0.7. It is 
the first analytical estimation that has obtained subdiffusion in the frozen 
turbulence. The problem is that this method can provide only the scaling of 
X, but not the statistics of the test particle trajectories. 

Important analytical results in the study of this special case were obtained 
in the last decades by developing new statistical methods (the decorrelation 
traj ectory method DTM |15| and the nested subensemble approach NSA 
fl6[). They determine the LVC, the time dependent diffusion coefficients 
Di(t), the probability of displacements until decorrelation and other statisti¬ 
cal averages. It was shown that the presence of trapping determines memory 
effects in L(t) and a rich class of anomalous diffusion regimes in the presence 


of a decorrelation mechanism 16]. The trapping has also collective effects. 


It determines coherence in the stochastic motion in the sense that bundles 
of neighboring trajectories form localized structures similar to fluid vortices. 

DTM reduces the problem of determining the statistical behavior of the 
stochastic trajectories to the calculation of weighted averages of smooth, 
deterministic trajectories obtained from the Eulerian correlation of the po¬ 


tential 15]. DTM is a semi-analytical statistical approach that satisfies the 


statistical consequences of the invariance of the potential. 
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The main idea of our approach is to study the stochastic equations (P 
in subensembles of realizations of the stochastic field, which contain all re¬ 
alizations that have the same values of the stochastic potential and velocity 
in the starting point of the trajectories 

(S): 0(0,0)=/, v (0,0) = v°. (11) 


A average potential <f> 5 and velocity V 5 exist in each subensemble. They are 
space-time dependent functions determined by the EC of the potential as 


0 £(M) 0 E y (x,t) _ 0 W) 

’ £7(o,o) x v 2 y v y 2 ’ 


( 12 ) 


V s (x,l)=e 2 xVl s (x,t) (13) 

where x = ( x , y ), 0°, v f are the parameters of the subensemble and the sub¬ 
scripts represent derivatives (P) = dE/dxi). The amplitudes of the velocities 
are V% — —E yy (0), Vy — —E xx (0). 

The stochastic equations are studied in each subensemble (S), where tra¬ 
jectories with a high degree of similarity are obtained due to the supplemen¬ 
tary initial conditions (TlTT) . Neglecting the fluctuations of the trajectories, 
the average trajectory in (S) (the decorrelation trajectory) is obtained by 
averaging Eqs. (0Q) in (S). One obtains average equations with the same 
structure as the equations in each realization (|T|) , but having the stochastic 
potential replaced by the subensemble (conditional) average potential 


dX 5 

dt 


- e z x V$ + V^e,,. 


(14) 


This approximation is validated in 16], where it is shown that the DTM is 


the first order in a systematic expansion, and that the results obtained in 
the second order are close to those of the first order. 

The Lagrangian correlations are obtained as averages over the decorre¬ 
lation trajectories, by summing the contribution of each subensemble (S), 
weighted by the probability that a realization belongs to the subensemble. 
In particular, the LVC is approximated using the decorrelation trajectories 


as 


L u (t) ^ I d0 u P(0 u ) / d/du 2 °P(/)PK) v. 


, 0 n „,0 


dXf(t) 

dt 


(15) 









where P is the Gaussian probability 


P(0°) 



om 

2/^ 2 / ’ 


p K°) 



(c°n 

2^ 2 J • 


Thus, the DTM is essentially an analytical method, which determines the 
diffusion coefficients from a set of simple, deterministic trajectories that result 
from the EC of the stochastic held. However, in most cases, the DT’s have 
to be numerically integrated. A computer code was developed for calculating 
the time dependent diffusion coefficients. It determines the DT’s (TO) for a 
large enough number of subensembles and performs the integrals in Eq. (11511 . 
These computer calculations are at PC level, for times of the order of 10 
minutes. 


3. Effects of trajectory trapping on test particle statistics 


3.1. Memory effects and anomalous diffusion coefficients 

The Lagrangian velocity correlation (LVC) determines the time depen¬ 
dent diffusion coefficient (jHj) and the mean square displacement ([TJ) . It is also 
a measure of the statistical memory of the stochastic motion. 

The presence of trajectory trapping leads to a specific shape of the LVC 
for a frozen turbulence 0(x). This function decays to zero in a time of the 
order tji but at later times it becomes negative, it reaches a minimum and 
then it decays to zero having a long, negative tail. The tail has a power law 
decay with an exponent that depends on the EC of the potential [35J. The 
positive and negative parts compensate such that the integral of L(t), the 
running diffusion coefficient D(t), decays to zero. The transport in such two- 
dimensional potential is thus subdiffusive. The tail of the LVC shows that the 
stochastic trajectories in frozen turbulence have long time memory. Thus, the 
LVC for incompressible two-dimensional velocity fields is completely different 
from the LVC of diffusion without trapping. The latter is decaying from V 2 
to zero in a time of the order r //. 

This stochastic process is unstable in the sense that any weak perturba¬ 
tion produces a strong influence on the transport. A perturbation represents 
a decorrelation mechanism and its strength is characterized by a decorre¬ 
lation time Tfi- The weak perturbations produce long decorrelation times, 
Td 3> Tfi. In the absence of trapping, such a weak perturbation does not 
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produce a modification of the diffusion coefficient because the LVC is zero at 
t > Tfi In the presence of trapping, which is characterized by long time LVC, 
such perturbation influences the tail of the LVC and destroys the equilibrium 
between the positive and the negative parts. Consequently, the diffusion co¬ 
efficient is a decreasing function of r d . It means that when the decorrelation 
mechanism becomes weaker {r d increases) the transport decreases. This is a 
consequence of the fact that the long time LVC is negative. This behavior 
is completely different from that obtained in stochastic fields that do not 
produce trapping. In this case, the transport is stable to the weak perturba¬ 
tions, because the LVC is zero for t > Tfi. An influence of the decorrelation 
can appear only when the perturbation is strong such that r d < Tfi and it 
determines the increase of the diffusion coefficient with the increase of r d . 

The decorrelation can be produced for instance by the time variation of 
the stochastic potential, which destroys the Lagrangian correlations at t > r c . 
The transport becomes diffusive with an asymptotic diffusion coefficient that 
has the trapping scaling ()I0|1 . and thus it is a decreasing function of r c . This 
behavior is determined by the fact that a stronger perturbation (with smaller 
r d ) liberates a larger number of trapped trajectories, which contribute to the 
diffusion. For other types of perturbations, the interaction with the trapping 
process produces more complicated nonlinear effects. For instance, particle 
collisions lead to the generation of a positive bump on the tail of the LVC 


19] due to the property of the 2-dimensional Brownian motion of returning 


in the already visited places. 

The average velocity V d e y does not provide a decorrelation mechanism, 
but it determines an average potential that changes the structure of the 
contour lines by producing bunches of opened lines (as discussed in Section 
2.1). A strong modification of the diffusion coefficient is produced by the 
average velocity in the presence of trajectory trapping (A* > 1 and K > 1), 
which consists of a large increase of the parallel diffusion D y and of the 
decrease of the perpendicular diffusion D x [36], Also, an interesting process 
of ’’acceleration” appears. Since a fraction of the trajectories 7i tr are trapped, 
they cannot participate to the average flux, which is produced only by the 
free trajectories. The velocity field has zero divergence, thus the distribution 
of the Lagrangian velocity is time independent, and in particular, the average 
Lagrangian velocity is equal to V d . Due to this, the average velocity of the free 
particles Vf has to be lager that V d such that njVf = V d . Here rif — 1 — n tr 
is the fraction of free trajectories. The increased velocity Vj is determined 
by the selection of the stochastic velocity along the opened contour lines of 
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Figure 2: The probability of displacements along the average velocity V,ie y obtained with 
the DTM for K = 10 and r = Vd/V = 0.2 at several time moments. 


the potential, which are situated preferentially in regions where the Eulerian 
velocity is oriented along V^e^. The trapping determines in these conditions 
a strong modification of the distribution of displacements, which elongates 
in the direction of Vde y and eventually it splits in two parts (see Figure 2). 

The case of a stochastic potential that moves with the velocity V r ie y is 
obtained by a Galilean transformation. This leads to opposite flows: trapped 
particles move with the potential and the free ones move in the opposite 
direction such that the total flux is zero, n t rVd + rifVf = 0. 


3.2. Coherence and trajectory structures 

The statistical characteristics of the trapped and of the free trajectories 
were studied in 16]. We have shown that the two types of trajectories have 


completely different statistical characteristics. 

The trapped trajectories have a quasi-coherent behavior. Their aver¬ 
age displacement, dispersion and probability of displacements saturate in a 
time t s . The time evolution of the square distance between two trajectories 
( Sx 2 (t )) is very slow showing that neighboring particles have a coherent mo¬ 
tion for a long time, much longer than r s . They are characterized by a strong 
clump effect with the increase of ( 5x 2 {t )) that is slower than the Richardson 
law. These trajectories form quasi-coherent structures which are similar to 
fluid vortices and represent eddying regions. 

In time dependent potentials (with finite r c ), the structures with t s > r c 
are destroyed and the corresponding trajectories contribute to the diffusion 
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1.5 



Figure 3: The formation of quasi-coherent trajectory structures in a time independent 
stream function ( K,t c = oo). The time evolution of nt r , S x , S y for Vd = 0 (dashed lines) 
and Vd = 0.2 (continuous lines). 


process. These free trajectories have a continuously growing average displace¬ 
ment and dispersion. They have incoherent behavior and the clump effect 
is absent. The probability of short time displacements are non-Gaussian for 
both types of trajectories. 

The study of vorticity separation (presented in Section 4.1) requires more 
information on the trajectory structures than obtained in [16] . It is necessary 
to determine the fraction of trapped trajectories n tr and the average size of 
the structures Si as functions of the average velocity V0 

These statistical quantities are obtained using the DTM as weighted av¬ 
erages of the decorrelation trajectories X 5 (f) in the static potential. The 
solutions of Eq. (fT4]l are periodic functions in this case with the periods 
T(0°,v°) that depend on the subensemble. The fraction of trapped trajec¬ 
tories at time t is 

n tr (t)= j d0°P(0°) J dv°dv%P(v°)P(vi) c tr (t;(f) 0 , v°), (16) 

where Q r (t; 0°, v°) = 1 if t > T(0°, v°) and c tr {t ; 0°, v°) = 0 if t < T(0°, v°). 
The sizes of the trajectory structures on the i = x,y directions are 

Si{t) = f d0°P(0°) J dv°dv°P{vl)P(v°) A7 ax (t;0°,v°), (17) 
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Figure 4: The dependence of the parameters of the quasi-coherent structures on the average 
velocity. 

where 0°, v°) = rnax(Xf) if t > T(0°,v°) and 0°, v°) = 

0 if t < T(</>°,v°). These functions describe the growth of the trajectory 
structures. 

The time evolution of the fraction of trapped trajectories nt r {t ) and the 
average size of the structures S x (t), S y (t ) are plotted in Figure 3 for Vd — 0 
(dashed lines) and for Vd = 0.2 (continuous lines). For Vd = 0, the structures 
continuously grow due to the large size stream lines that lead to large periods 
of the corresponding trajectories. An average velocity generates bunches of 
opened lines of the total potential and limits the size of the islands of closed 
contour lines. The size of the structures and n tr {t) saturate at finite Vd, as 
seen in Figure 3. The average velocity determines the decrease of n tr and of 
the size of the structures in the perpendicular direction S x , while the parallel 
size S y is increased. The structures are elongated by the velocity. 

The saturation values of the parameters of the structures n tr , S x , S y 
and nf — 1 — n tr are plotted in Figure 4 as functions of Vd- One can see 
that the trajectory structures exists only for small average velocities with 
r = Vd/V < 1 (or AT > 1). When r > 1, the structures are destroyed by the 
average velocity, which opens all the contour lines of the total potential. 

4. Effects of trapping on turbulence evolution 

The idea of the strong relation between tracer trapping or eddying and 
turbulence evolution is sustained by analyzing two completely different phys- 
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ical systems: the relaxation of turbulence in two-dimensional ideal fluids and 
the drift turbulence in magnetically confined plasmas. Fluids are stable and 
turbulence is externally produced by steering, while inhomogeneous plasmas 
are unstable and generate spontaneously turbulence. However, both sys¬ 
tems are characterized by the inverse cascade of the energy in the evolution 
of turbulence and by the tendency to self-organization. We show that in 
both cases these processes are determined by trajectory trapping, but with 
different mechanisms, which are essentially determined by the different non- 
linearities of the systems. The case of fluid turbulence is analyzed in Section 
4.1, where the statistics of the trajectories of the vorticity elements in the 
neighborhood of a large vortex is studied. Section 4.2 contains a short review 
of recent results on the evolution of the drift turbulence in magnetically con¬ 
fined plasmas. It is introduced in order to underline the idea that trajectory 
trapping has the main role in both cases, despite the fact that the physical 
processes that determine the inverse cascade are completely different. 

4-1- Ideal fluids 

Ideal fluids are described by the Euler equation 

d t uj + v • Vu; = 0, (18) 

v = e z x V0, (19) 

u = A 0, (20) 

where uj is the vorticity, 0 is the stream function and v is the fluid velocity. 
Thus vorticity elements are advected by the velocity held and the vorticity is 
conserved on the trajectories. The nonlinearity of the process is determined 
by the relations (jT9|) - (l20]) between v and u, which show that vorticity is an 
active held. 

Numerical simulations have shown that the vorticity held evolves to orga¬ 
nized states characterized by isolated vorticity peaks with large amplitudes. 
A process of separation of positive and negative vorticity and vorticity con¬ 
centration occurs during the relaxation of turbulent states. Our aim is to 
understand how is possible to appear vorticity separation in a turbulent state. 

In order to study this tendency, we consider a turbulent state as initial 
condition of Eq. (USD that is represented by a stochastic stream function 
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with Gaussian distribution and with the Eulerian correlation A(x). We have 
taken for the results presented in the figures 



( 21 ) 


where r = |x|/A. This function is normalized to E( 0) = 1 and its space 
integral is zero. The latter condition is due to the choice of the total vorticity 
equal to zero. The presence of a large scale vortex is modeled by an average 
stream function (0(x)) = xV ( i that accounts for an average velocity Vde y . 
We consider first an uniform velocity, and the effect of a small transverse 
gradient of Vde y is discussed at the end of this section. 

Trajectory trapping occurs when A* > 1 and K > 1, as discussed in 
Section 3.1. The first condition implies large amplitudes of the stochastic 
component of the velocity compared to the average velocity V > Vd- The sec¬ 
ond condition corresponds to large correlation times of the stochastic stream 
function r c > r // = X c /V. The temporal evolution of the stream function 
is determined in ideal fluids by the motion of the vorticity elements, which 
modifies the vorticity field. This process is slow compared to the time of 
flight, which means that fluid turbulence has K > 1. Thus, trapping occurs 
in two-dimensional fluid turbulence when the average velocity is zero or small 


(Vd < V). 


The stream function and the vorticity are correlated and their two-point 
correlation function results from A(x) using Eq. (12Uj) 


E<j,u{x) = (</>(x') u;(x' + x)) = AA(x) 


( 22 ) 


Thus, the stream function and the vorticity in the same point (x = 0) have 
opposite signs since E( 0) is the maximum of A(x) and AA(0) < 0. This 
correlation, determined by the nonlinearity of the fluid flow, has an important 
role in vorticity separation. 

The trajectories of the vorticity elements are solutions of Eq. ©• Since 
they are independent on the vorticity and on its sign, it is not evident how 
vorticity separation can occur. 

We show that the average velocity Vd.e y determines two transverse flows 
with opposite directions. They appear in the presence of eddying and their 
direction is correlated with the sign of the stream function corresponding to 
the stream line on which the trajectory is trapped. 
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The average displacement of the vorticity elements is determined using 
DTM by summing the contributions of all subensembles. An equation similar 
to (fT5l) is obtained 

<x(f)> = I d(fp P(0°) J dv° P(v°) X 5 (f), (23) 

where X 5 (f) is the decorrelation trajectory in the subensemble S determined 
from Eq. (TT4l) . The component perpendicular on the average velocity is zero. 
Actually it is of the order 10 -5 , which is the accuracy of the numerical cal¬ 
culation. This is the correct result for an incompressible velocity field. The 
average Lagrangian velocity has to be equal to the Eulerian average veloc¬ 
ity, which has no component in the e x direction. However, the conditional 
average displacement for a given value of the initial stream function 0° 


(!(()),. = P(<p°) / <iv° P(v°) X s (t) 


(24) 


is not zero, but it is an anti-symmetrical function of 0°, as shown in Figure 5. 
One can see that the displacements are correlated with the stream function 
and that (x( t ))^ 0 and of 0° have the same sign. The displacements are 
negligible at small time when there is no trapping. As time increases, the 


maximum of 


'At)), 


increases and eventually saturates. 

Due to the anti-correlation of the stream function and vorticity in the 
same point (£0,^(0) < 0), the positive peaks of the stream function corre¬ 
spond to negative vorticity, while the negative ones correspond to positive 
vorticity. Thus the correlation of the perpendicular displacement with the 
stream function corresponds to a correlation between displacements and vor¬ 
ticity. 

The total average displacements for positive 0°, ( x(t)) + , and respectively 
for negative 0°, ( x(t))_ , are represented in Figure 6. They are symmetrical 
and compensate, as required by the zero divergence of the velocity field. 
These average displacements determine opposite average velocities in the 
case of time dependent vorticity fields 


^+ = 


Ar c )). 


T r 


(x(t c )) 

Vs- = = -V,+, 


T r 


(25) 


where r r is the correlation time. 


16 






15 



-15 - 1 - 1 - 1 - 1 - 1 - 

- 3 - 2-1 0 1 2 3 


Figure 5: The average perpendicular displacement as function on the initial stream func¬ 
tion for r = Vd/V = 0.2 and K = oo. 



Figure 6: The average perpendicular displacements for the positive and negative vorticity 
(r = Vd/V = 0.2 and K = oo). 


17 









0.12 



Figure 7: The average normalized velocity of positive and negative vorticity separation, 
V a /V as function of the Kubo number for r = 0.2. 


In terms of vorticity, the negative vorticity elements that are trapped on 
the stream lines with 0° > 0 move with the velocity V s+ in the direction 
of the gradient of the large scale stream function. The positive vorticity 
elements have an opposite average velocity. Thus, the positive and negative 
vorticity elements that evolve on trapped trajectories separate with a speed 
V s = V s+ — V s - = 2V S+ . This speed is plotted in Figure 7 as function of 
the Kubo number K = r c /ryy One can see that V s is maximum for K of 
the order one, and it decays at K > 1 as K due to the saturation of 
the average displacements. The maximum value is smaller than 10, but not 
negligible. 

The physical mechanism for the vorticity separation consists of the influ¬ 
ence of the average velocity on the small structures of trapped trajectories. 
Vd determines a difference between the average velocities on the two sides of 
the closed trajectories that are oriented along V0e y . The velocity is d x (j) + V0 
on one side and d x (j) — V ( j on the other side. The difference of velocities leads 
to the accumulation of the vorticity elements on the small velocity side and to 
the increase of the vorticity in that region. As the direction of the velocity on 
the contours of 0 depend on the sign of 0, the accumulation of the vorticity 
elements occurs on opposite sides for contours corresponding to 0 and —0. 
The time variation of the velocity held destroys the trajectory structures in a 
time of the order r c , and the coherent displacements of the vorticity elements 
produce transverse average velocities. 


18 




The average displacement (y(t)) along the average velocity Vd& y is pro¬ 
duced, as explained in Section 3.1, only by the free trajectories. They de¬ 
termine the elongation and eventually the splitting of the distribution of 
displacements (see Figure 2). Thus, the vorticity elements spread along Vde y 
while they have transverse drifts. 

The average flows of the vorticity elements determine the modification 
of the vorticity distribution. Due to the correlation (J22]), this changes the 
distribution of the stream function: the contour lines of the stream function 
follow the average motion of the vorticity. 

These findings are in agreement with the image provided by the numerical 
simulations (see for instance 37]), which show the emergence of large vortices 
that elongate and attract small vortices of the same sign. We note that the 
process analyzed here is different from the vortex drift transverse to a sheared 
flow analyzed in [38] since it involves vorticity elements in turbulent flows 
and not localized vortices. In our case the drift is produced by an average 
velocity and the existence of a gradient of the background vorticity is not 
necessary. 

A vortex determines a velocity that has a gradient in the transverse di¬ 
rection. The gradient length is much larger than the correlation length of 
the stream function in the case of a large vortex in a small scale stochastic 
velocity held. In these conditions, the effects of the gradient can be esti¬ 
mated using the dependence on Vd of the statistical characteristics obtained 
at constant Vd- 

The transverse drift of the small structures with the same sign of the 
vorticity as the large vortex is in the direction of the gradient and brings the 
structure toward larger Vd- This determines the decrease of the transverse 
dimension of the structure and the increase of its elongation, as seen in 
Figure 4. When Vd > V, the small structures are destroyed. The fraction 


of trapped trajectories and the size of the structures decrease rapidly for 
Vd > V Their vorticity is included in the large vortex, and it leads to the 
increase of the large scale energy (inverse cascade). On the contrary, the 
structures with opposite sign of vorticity drift toward smaller Vd, against the 
gradient. Their size increases and their elongation is reduced. 

In conclusion, a large size vortex in a decaying two-dimensional turbu¬ 
lence attracts the small size vortices with the same sign of vorticity while the 
small vortices of the opposite sign are repelled. This effect of vorticity sepa¬ 
ration according to its sign is produced by trapping of the vorticity elements 
combined with the existence of the correlation of the stream function with 
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the vorticity imposed by the nonlinearity (1201) . 


4-2. Drift turbulence 

The nonlinear stage in the evolution of drift type turbulence is essentially 
a consequence of ion trajectory trapping or eddying in the structure of the 
stochastic potential. This conclusion is drawn from a study of test modes 
on turbulent plasmas [2f 


39] based on a new Lagrangian approach that 


extends the type of methods initiated by Dupree 40] to the nonlinear regime 
characterized by trapping. Drift waves are unstable due to the electron ki¬ 
netic effects that produce the dissipation mechanism to release the energy, 
combined with the ion polarization drift |41]. Beside this, the polarization 
drift has a more complex influence determined by its nonzero divergence, 
which produces a weak compressibility effect in the background turbulence. 

Test modes on turbulent plasmas were studied for drift turbulence in con¬ 
stant magnetic field starting from the basic description provided by the drift 
kinetic equations. Analytical expressions are derived, which approximate the 
growth rate y(k) and the frequency cj(k) of a test mode with wave number 
k = (k x , k y ) as functions of the characteristics of the background turbulence. 
These functions provide the tendency of turbulence evolution. 

The dispersion relation of the test modes in turbulent plasma is shown 
to be the same as in quiescent plasma, except for a time dependent func¬ 
tion M(t ), which embeds all the effects of the background turbulence. This 
function appears in the propagator of the test mode, which is obtained using 
the formal solution of the ion equation with the method of characteristics. 
The characteristics are ion trajectories in the stochastic background poten¬ 
tial. They are described by equations similar to Eq. ©. with the difference 
that the potential </> is moving with the effective diamagnetic velocity V^e^, 
which leads to the replacement of y with y' = y — V&t. The function M(t) is 
evaluated using the DTM, and the dispersion relation is analytically solved. 

Drift turbulence develops in the initial stage on a wide range of wave 
numbers. Ion trajectories are not trapped at these small amplitudes of the 
background turbulence and they have Gaussian distribution. Their diffusion 
determines the damping of the large wave number modes. The correlation 
lengths Xi remain close to the ion Larmor radius p, L during this stage. Tur¬ 
bulence amplitude /3 increases and the shape of the EC is not changed. 

When the amplitude reaches values that make AT > 1, ion trajectory 
trapping appears and generates vortical structures of trapped ions as dis¬ 
cussed in Section 3.2. They determine the decrease of the frequencies, which 
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leads to the decrease of 7 and to the displacement of the unstable range 
of wave numbers to smaller values. The maximum of the growing rate is 
displaced from k y = 1 /p i to k y = 1/S, where S is the average size of tra¬ 
jectory structures. Turbulence amplitude continues to increase in this stage, 
but with a smaller rate. This determines the increase of the fraction of the 
trapped ions n tr and of the size S of the trajectory structures. 

Large scale correlations are generated through a nonlinear process of de¬ 
creasing the wave number of the unstable modes. In the same time, the 
diffusive damping produced by the free trajectories continues to act on the 
large k modes, which are attenuated. This process of inverse cascade of 
turbulence energy stops well before the size of the potential cell becomes 
comparable with the system size. A different effect appears at larger trap¬ 
ping rates, which is also connected with trapping. 

When the fraction of trapped ions nt r becomes comparable with the frac¬ 
tion of free ions rif, ion flows are generated by the trapping process as ex¬ 
plained in Section 3.1. Trapped ions move with the potential, while free 
ions move in the opposite direction such that the total flux is zero. This 
determines the splitting of the distribution of ion displacement (as in Figure 
2 ). As shown in [25[, an essential change of test modes is produced. The 
drift modes are damped and a different type of modes, zonal flow modes, 
are generated in this strongly nonlinear regime. They have k y = 0 and very 
small frequencies and are produced by the combined action of the ion flows 
and of the compressibility due to the polarization drift in the background 
turbulence. The damping of the drift modes is not determined by the zonal 
flow modes as generally believed. There is only an indirect contribution 
through the diffusive damping, which is increased by the zonal flow modes 
|39j. The decay of the drift turbulence determines the decrease of n tr and of 
the growth rate of the zonal flow modes, which is proportional to n tr . Con¬ 
sequently, drift turbulence does not saturate, but it oscillates between weak 
and strong trapping. The predator-prey paradigm is not sustained by these 
results, although there is time correlation between the growth of zonal flow 
modes and the damping of the drift modes. 

In conclusion, trapping determines in the first phase the increase of the 
correlation length by displacing the range of unstable modes to small k ~ 
1/S. In the second stage, when it becomes stronger and produces ion flows in 
the moving background turbulence, trapping leads to the nonlinear damping 
of the drift modes. In the same time, trapping combined with the polarization 
drift generates zonal flow modes. 
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5. Conclusion 


Particle motion in two-dimensional incompressible stochastic velocity fields 
is characterized by a mixture of random and quasi-coherent aspects, repre¬ 
sented by a random sequence of large jumps and trapping or eddying events. 
Trajectory trapping generates quasi-coherent structures of trajectories. 

The development of statistical methods that are adequate for the descrip¬ 
tion of such complex trajectories has contributed to the understanding of the 
nonstandard characteristics of the transport in these systems. Various as¬ 
pects of the turbulent transport were analyzed in the test particle approach 
using these semi analytical methods. These studies are now in a new stage 
that consists of developing Lagrangian methods for the study of the evolu¬ 
tion of the turbulence in the nonlinear regime. First results are reported in 


25j and in the present paper. They lead to the conclusion that trajectory 
trapping has the essential role in turbulence evolution. The physical mech¬ 
anism depends on the specific nonlinearity of each system, but it is related 
to trajectory trapping. It appears that the quasi-coherent component of the 
trajectories determine the nonlinear effects in turbulence while the random 
component provides a damping or dissipation process. 

In the case of fluids, the trapped vorticity elements have an average 
motion, which leads to average velocities with the directions dependent on 
the vorticity sign. Then, due to the correlation E^, the stream function 
is also transported by the perpendicular average flows. The effect of this 
process is the separation of the vorticity with different signs since negative 
vorticity moves in the sense of the stream function gradient and positive 
vorticity in the opposite sense. A small gradient of the average velocity that 
characterizes a large scale vortex determines the decrease of the attracted 
small vortices and eventually the inclusion of their vorticity in the large 
vortex. This process determines the flow of the energy from small to large 
scales. 

In the case of drift turbulence, the trapped ions determine the displace¬ 
ment of the wave numbers of the most unstable modes towards small values, 
while the free ions contribute to the diffusive damping of the large k modes. 
The correlation length of the stochastic potential increases until the process 
of nonlinear damping becomes dominant. Ion trapping has an important role 
in this process and also in the generation of the zonal flow modes. 
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